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Using the hierarchy of scales between the mass, M, and the width, T, of a heavy, unstable particle 
we construct an effective theory that allows calculations for resonant processes to be systematically 
expanded in powers of the coupling a and T/M. We illustrate the method by computing the next- 
to-leading order line shape of a scalar resonance in an abelian gauge- Yukawa model. 
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Higher-order calculations for processes involving mas- 
sive, unstable particles close to resonance suffer from the 
breakdown of ordinary perturbation theory, since the in- 
termediate propagator becomes singular. This singular- 
ity is avoided if the finite width, T, of the unstable parti- 
cle is taken into account in the construction of the prop- 
agator via resummation of self-energy insertions. There 
are a number of approaches along this line to avoid the 
problem However, so far there is no method that 
allows to systematically improve the accuracy of calcula- 
tions order by order in perturbation theory. The purpose 
of this letter is to present such a method. 

We are concerned with processes involving an unstable 
particle close to resonance. The main idea is to exploit 
the hierarchy of scales T <C M, where M is the pole mass, 
in order to systematically organize the calculations in a 
series in the coupling, a, and T/M. While the expan- 
sion in a is standard, we construct an effective theory 
to perform the expansion in T/M. A first step in this 
direction has been presented in |2|. The main idea of 
our approach is similar to non-relativistic QCD, where 
an expansion in a and the velocity of the heavy quarks is 
made. We will identify all relevant modes and use them 
to write the operators of the Lagrangian of the effective 
theory. This Lagrangian is then matched to the underly- 
ing theory, using the method of regions Q . In this letter 
we will outline the basic idea and we refer to for more 
details. 

We consider a toy model that involves a massive scalar 
field, 4>, and two fermion fields. The scalar as well as 
one of the fermion fields, ip, (the "electron") are charged 
under an abelian gauge symmetry, whereas the other 
fermion, \i (the "neutrino") is neutral. The model al- 
lows for the scalar to decay into an electron-neutrino pair 
through a Yukawa interaction. The Lagrangian is 



C 



+ vH>x + y*<Pxi> - -^{4>\ 



(1) 



counterterm Lagrangian and D 



u g = 9 7(4tt), 



(1W*)/(4tt) 



d — igA. We define 
(at the scale fx) and 



assume a g ~ a y ~ a, and ot\ = X/(4tt) ~ o? /(JTk). 

We would like to obtain the totally inclusive cross sec- 
tion for the process 



9(q) + e (p)->X 



(2) 



as a function of s = (p + q) 2 by calculating the forward 
scattering amplitude T(s) and taking its imaginary part. 
(The total cross section of process (0 has an initial state 
collinear singularity which has to be absorbed into the 
electron distribution function. In what follows it is un- 
derstood that this singularity is subtracted minimally.) 
In particular, we are interested in the region s « M 2 , 
or more precisely s — M 2 ~ MT ~ aM 2 <C M 2 . In 
this kinematic region the cross section is enhanced due 
to the propagator of the scalar. Furthermore, at each 
order in a we get additional contributions proportional 
to aM 2 /(s — M 2 ) ~ 1 due to self-energy insertions. 

We now turn to the main part of this letter and dis- 
cuss how to construct the effective theory. Our approach 
is based on the hierarchy of scales T -C M . Thus, we 
systematically expand the cross section in powers of a 
and 



5 = 



s-M 2 
M 2 



M' 



(3) 



where M and £ c t denote the renormalized mass and the 



In a theory that formulates this expansion correctly, 
other issues like resummation of self-energy insertions 
and gauge invariance are taken care of automatically. 
In a first step we integrate out hard momenta k ~ M. 
The effective theory will then not contain any longer dy- 
namical hard modes since their effect is included in the 
coefficients of the operators. The hard effects are associ- 
ated with what is usually called factorizable corrections, 
whereas the effects of the dynamical modes correspond 
to the non- factorizable corrections 0- On the level of 
Feynman diagrams, this amounts to using the method of 
regions to separate loop integrals into various contribu- 
tions 0. The hard part is obtained by expanding the 
integrand in S. The difference between the full integral 



2 



and its hard part has to be reproduced by modes cor- 
responding to momentum configurations that are near 
mass-shell. The main task is to identify these modes, 
and to write the operators of the effective Lagrangian 
in terms of the corresponding field operators and then 
to compute the coefficients of the operators by matching 
(up to a certain order in a and 8). 

Our goal is to carry out this programme for our model 
to an order in a and 5 that is sufficient to compute 
7~(o) _|_ the forward scattering amplitude at next- 

to-leading order (NLO), where 7~(°) sums up all terms 
that scale as (a/S) n ~ 1 and 7"W contains all terms that 
are suppressed by an additional power of a or S. 

The basic process under consideration is the following: 
we start with highly energetic fermions, produce a near 
mass-shell scalar which then decays again into highly en- 
ergetic fermions. Accordingly we split the effective La- 
grangian into three parts. Roughly speaking, the first, 
£hset, describes the heavy scalar field near mass-shell 
and its interaction with the gauge field. The second part, 
£scet, describes energetic (charged) fermions and their 
interactions with the gauge field. Finally, the third part, 
£; n t , describes the external fermions and how they inter- 
act to produce the final state. We will discuss these three 
parts in turn. 

The construction of £hset follows closely the construc- 
tion of the effective Lagrangian for heavy quark effec- 
tive theory (HQET) [jj. We write the momentum of the 
scalar particle near resonance as P = Mv + k, where the 
velocity vector v satisfies v 2 = 1 and the residual mo- 
mentum k scales as MS. We will call such a scalar field 
a "soft" field (in the term "resonant" has been used). 
Thus, for a soft scalar field we have P 2 — M 2 ~ MS and 
this remains true if the scalar particle interacts with a soft 
gauge boson with momentum MS. In analogy to HQET 
we remove the rapid spatial variation e ~ tMv ' x from the 
<A-ficld and define 



scales as vk ~ M. For the other we find 



(j) v (x) 



AMv-x 



(4) 



where P+ projects onto the positive frequency part to 
ensure that <p v is a pure destruction field. We now write 
the effective Lagrangian in terms of <p v and construct the 
bilinear terms so as to reproduce the two-point function 
close to resonance. Denoting the complex pole of the 
propagator by s and the residue at the pole by the 
propagator can be written as 



i R4 
P 2- 



iR d 



2Mvk + k 2 - (s - M 2 ) 



We define the matching coefficient 

s - M 2 



M 



(5) 



(6) 



and Oj = - 
solutions to P 2 



(va) for any vector. There are two 
= s, one of which is irrelevant since it 



vk = -M + yjM 2 + MA - k 2 
A_tf_ 

7 



4k 2 



8M 



0(S 3 ), 



(7) 



where we expanded in S in the second line. Therefore, 
the bilinear terms are given by 



C-rhch — 



2M4 [iv-D s - = 



2McPl 



2M 



A 2 \ 
8MJ 



(8) 



where D s = d — igA s denotes the soft covariant deriva- 
tive. In obtaining we exploited the fact that the 
gauge invariance of the full Lagrangian is not broken by 
the separation into hard and soft parts. Therefore, the ef- 
fective Lagrangian must be gauge invariant as well and we 
can obtain the interaction of the scalar with the soft pho- 
ton simply by replacing d — > D s . The gauge invariance of 
A follows from the gauge invariance of s and M. Further- 
more, A is given entirely by hard contributions, which 
justifies its interpretation as matching coefficient. Using 
© we can express it in terms of the hard part of the self- 
energy Tl h (s). Writing IL h (s) = M 2 E^^lL^), where 
it is understood that n^'^ ~ a k , we obtain 



A = ^A« = 

i 

tfnM + AfnM+nfWnti 



(9) 



0) 



Explicit results for A^ and A^ 2 ) in the MS and pole 
renormalization scheme can be found in Here we 
only note that in the pole scheme s = M 2 — iMT, so 
A = — iT when M = M. Inserting the expansion © 
into JSJl and supplementing with the kinetic terms 
for soft photons and fermions we obtain 



-HSET 



2M0t Uv ■ D s - —j cf,. 



2M<pl 



f { (iDs.t) 2 | [A' 1 '] 2 A^ 



2M 



8M 



(10) 



Each term in £hset can be assigned a scaling power in S. 
In momentum space the propagator of the <p v field scales 
as 1/8. Hence, because / d 4 k counts as 8 4 , the soft scalar 
field <t> v (x) scales as S 3 / 2 . Since A^ - D s ~ MS, both 
terms in the first line of ((HJ) scale as 5 A and are leading 
terms. The terms in the second line are suppressed by 
one power in S or a. Finally, since j4£scales as 5 and the 
soft fermion fields scale as <5 3 / 2 (see |7|) the terms in the 
last line of (jHJ scale as S 4 . In JSJ we have left out terms 



3 



further suppressed in 5 or a. As we will see, they are 
not needed for the calculation of the line shape at NLO. 
However, we stress that the expansion can be performed 
to whatever accuracy is needed. 

We note that computing the scalar propagator to all 
orders in S using £hset does not reproduce iJSJ. Instead 
near resonance we obtain ivu^ 1 R C S4> / '{P 2 — s), where 
uj- 1 = {M 2 + MA-k\) l ' 2 /M = 1 + O(d,o). The differ- 
ence in the normalization is taken into account in match- 
ing calculations by an additional wave-function normal- 
ization factor nj -1 / 2 for each external cf> v -\me in the ef- 
fective theory. 

Next, we turn to the construction of the effective La- 
grangian, £scet, associated with the energetic fermions. 
We need a "collinear" mode to describe a fermion with 
large momentum in the say ft— direction. Such modes 
have been discussed previously within the context of soft- 
collinear effective theory (SCET) (a. The Lagrangian 
has been worked out to order 8 in Q and we can take 
the parts relevant to us from there. (What we call "soft" 
here what is usually called "ultrasoft" in the context of 
SCET and in the power counting our 5 corresponds to A 2 
in 0.) For each direction defined by an energetic par- 
ticle we introduce two reference light-like vectors, n±, 
with n 2 ^ = n 2 _ = and n + rt_ = 2 and we write the 
corresponding momentum as 

P" = (n+?)-f +j£ + (n-P)-p (11) 

where n + p ~ M, n^p ~ MS and p± ~ MS 1 ! 2 . Given a 
certain direction n_ we introduce the collinear field tp c 
which satisfies fi-ip c — 0. The terms relevant for the 
calculation of T*- -* + are then given by 

£scet = tpc (in-D + ipA_c- * , . iP±c J 4r V'c ■ 

\ in + U c + le ) 2 

(12) 

Since we are concerned with the forward scattering am- 
plitude, the only directions defined by energetic particles 
are given by the incoming electron and (anti)ncutrino. 
Thus, we have two sets of collinear modes, one for 
the incoming electron, ip c i, an d one for the incoming 
(anti)ncutrino, Xc2- Of course, in the case of the neutrino, 
the covariant derivatives in Ijl2|) have to be replaced by 
ordinary derivatives. All terms in (|12fl scale as S 2 . Terms 
of order S 5 ^ 2 and S 3 exist, but they are not needed for 
our application, since they would result in contributions 
suppressed by an additional power of a and, therefore, 
contribute only at NNLO. Again there is no difficulty in 
going to higher orders in the expansion if needed. 

The last part to consider is C ln t. It has to include 
operators that allow the production and decay of the un- 
stable particle. Without introducing additional modes it 
is not possible to include such vertices as ordinary inter- 
action terms in the effective Lagrangian Q . The reason is 
that the momenta associated with generic collinear fields 



ipd and \c2 do not add up to a momentum of the form 
P = Mv + k. Either we have to implement this kinematic 
constraint on our external states by hand 4] or we have 
to introduce a new "external-collinear" mode. Adopting 
the second option, we define an external-collinear mode 
with large momentum in the ft- direction by assigning it 
a momentum Mn_/2 + k, where k ~ S. This mode has 
the same virtuality MS 1 / 2 as a generic collinear mode but 
the momentum is not given by (fTTjl , because it has a fixed 
large component such that the two incoming fermions 
produce a scalar near mass shell. For such a mode it is 
useful to extract the fixed large momentum and to define 

1> n _(x) = e^^ n -^V + i> cl {x), (13) 

and similarly for Xn + ■ For the purpose of computing 
7~(o) -i-Xt 1 ) it is sufficient to take the first term of £scet, 
(|12fl . with a soft photon only to describe the interaction 
of the external-collinear fermions with the photons 

£± = ipnJn-D s —ip n __ +Xn+in + d — Xn+- (14) 

With the external-collinear modes we can implement 
the production and decay vertices as interaction terms 
in £i n t. It is also convenient to integrate out generic 
collinear fields and keep only the external-collinear modes 
in the effective theory. Because adding soft fields results 
in a further suppression in S we then find that we can 
restrict ourselves to 

Ant = Cy (f) v ^ n _ Xn+ + C y*4>tXn+^n- 

+ F^(^„_X„ + )(Xn + ^n_), (15) 

where C = l + 0(a) and F are the matching coefficients. 
The external fields scale as J 3 / 2 . Thus, an insertion of 
a <fiipx operator results in J d^x (f>vipn—Xn+ ^ 

5 1 ' 2 . The 

forward scattering amplitude can be obtained by two in- 
sertions of this operator. Taking into account the scaling 
of the external state (Pe~| ~ we see that T(°) - a/S. 
The four-fermion operator is suppressed in S and results 
in a contribution of order a to T. Thus, to compute T^ 1 ' 
we need the O(a) contribution to the matching co- 
efficient C, while F is only needed at tree level. 

The coefficient is obtained by matching the on- 
shell three-point function of a scalar field, an electron and 
a neutrino at order ya and at leading order in S. In par- 
ticular, this involves the computation of (the hard part) 
of the vertex diagram, and the additional wave-function 
normalization factor V3~ x l 2 mentioned above has to be 
taken into account. For the precise matching equation as 
well as the explicit result for we refer to Q. Here it 
suffices to say that these are standard loop calculations. 
To obtain F^ (the LO contribution to F) we have to 
match the four-point function at tree level, but include 
subleading terms in 5. The explicit result is F^ = 1/4. 
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FIG. 1: The line shape (in GeV" 2 ) in the effective the- 
ory at LO (light grey/magenta dashed) and NLO (light 
grey/magenta) and the LO cross section off resonance in the 
full theory (dark grey/blue dashed) as a function of the center- 
of-mass energy (in GeV). 



We have now completed the construction of the effec- 
tive Lagrangian C g = £hset +£± + A n t to an accuracy 
sufficient to compute T at NLO. At leading order there 
is only one diagram, involving two three-point vertices 
and one resonant scalar propagator. We get 



iT (0) = 



-iyy 

2MV 



[u(p)v(q)] [v(q)u(p)], 



(16) 



where we defined V = y/s — M — AW/2. In the effective 
theory there are three classes of diagrams that contribute 
to T^ 1 ) . Firstly, there are hard corrections consisting of a 
propagator insertion [AW] 2 /4 — AIA^ as well as a ver- 
tex insertion Secondly, there is a four-point vertex 
diagram due to the (ipx) (xi>) operator in £- lnt . The third 
class are soft-photon loop diagrams, corresponding to the 
non-factorizable corrections. Adding up all these contri- 
butions and using the explicit result for C^ 1 ' (in the MS 
scheme) we obtain 



; T (i 



: T (0). 



(17) 



31n- 



-2MV 



4 hi- 



-2MV 
M 2 



■ In ■ 



-2MV 



7 In 



2 In 



-2MV 




8T>M 



V 

2M 



where dj = on I (47r) . The initial state collinear singular- 
ities have been subtracted minimally and we denote the 
corresponding factorization scale by v to distinguish it 
from the renormalization scale fi. 

We can now perform the polarization average and take 
the imaginary part of (T(°) + TW)/s. This result de- 
scribes the line shape near resonance with a relative er- 
ror of a 2 . Moving away from the resonance, the relative 



error becomes of order unity, since 6 is not small any 
longer. To obtain a good description for all values of y/s, 
the result of the effective theory has to be matched to 
the off-resonance result of the full theory. 

In Figure ^ we show the leading order line shape in 
the effective theory and the tree-level (order a 2 ) cross 
section off resonance in the full theory. The two results 
agree in an intermediate region where both calculations 
are valid. This allows to obtain a consistent LO result 
for all values of y/s. We also show the NLO line shape. 
For the numerical results we have chosen to use the MS 
scheme with a y — a g — 0.1 and a>\ = (0.1) 2 /(47r). The 
pole mass is assumed to be M — 100 GeV which results 
in the MS value M = 98.8 GeV for the LO result and 
M = 99.1 GeV for the NLO result. Furthermore, we have 
chosen a variable factorization scale such that there are 
no large logarithms involving v. We remark that in order 
to obtain an improved NLO result for the whole region 
of y/s, the NLO line shape would have to be matched to 
the NLO off-resonance cross section in the full theory. 

The example considered here is based on a rather sim- 
ple toy model. Nevertheless, it allows to address the 
conceptual issues related to unstable particles. The main 
result is that, using an effective theory approach, calcula- 
tions can be performed in a systematic way in expanding 
in the small quantities a and T/M. Applying our method 
to the Standard Model might require more tedious calcu- 
lations, but the main result remains valid. In particular, 
as discussed in NNLO line-shape calculations now 
appear feasible. 
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